The role of structural evolution on the quantum conductance behavior of gold 

nanowires during stretching 
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Gold nanowires generated by mechanical stretching have been shown to adopt only three kinds of 
configurations where their atomic arrangements adjust such that either the [100], [111] or [110] zone 
axes lie parallel to the elongation direction. We have analyzed the relationship between structural 
rearrangements and electronic transport behavior during the elongation of Au nanowires for each 
of the three possibilities. We have used two independent experiments to tackle this problem, high 
resolution transmission high resolution electron microscopy to observe the atomic structure and a 
mechanically controlled break junction to measure the transport properties. We have estimated 
the conductance of nanowires using a theoretical method based on the extended Hiickel theory 
that takes into account the atom species and their positions. Aided by these calculations, we have 
consistently connected both sets of experimental results and modeled the evolution process of gold 
nanowires whose conductance lies within the first and third conductance quanta. We have also 
presented evidence that carbon acts as a contaminant, lowering the conductance of one-atom-thick 



PACS numbers: 68.37.Lp, 73.63.-b, 72.10.-d 



-a 

G 

o 
o 



> 

(N 

O 
<N 
O 
-i— > 



i 

G 

O 

o 



X 



I. INTRODUCTION 

The electron transport through nanometric conduc- 
tors attract a huge interest due to constant shrinkage of 
microelectronic devices.EJ In particular, metal nanowires 
(NW) display interesting quantum conductance behav- 
ior even at room temperature.! 2 ] From a practical point of 
view, NW's can be easily generated by putting in contact 
two metal surfaces, which are subsequently pulled apart. 
During the NW elongation and just before rupture, the 
conductance displays flat plateaus and abrupt jumps, 
which for metals such as Au, take a value of approxi- 
mately one conductance quantum Go = 2e 2 /h (where e 
is the electron charge and h is Planck's constant). 

In spite of the simplicity of the experimental pro- 
cedure, a new structure with a different evolution is 
observed for each NW generation and all conductance 
curves have plateaus and jumps, but they display dis- 
parate profiles.0 In order to overcome this difficulty, a 
simple statistical method has been usually applied to an- 
alyze the average behavior. Instead of considering the 
conductance as a function of the elongation, the trans- 
port properties can be represented as a histogram of con- 
ductance occurrence, in such a way that a flat plateau 
generates a histogram peakE By linearly adding the his- 
tograms associated to each conductance curve, a global 
histogram is generated, which describes the general ten- 
dencies of an ensemble of NW's. The global histogram 
displays well defined peaks close to the integer multiples 
of the conductance quantum; this fact has been adopted 
as a proof of the tendency to conductance quantization 
in metal NW'sH 

The statistical method, discussed above, provides in- 
formation on the average behavior but it hinders the 
study of NW transport properties. For example, it is not 



possible to get the detailed information on how structural 
factors influence the conductance evolution. For several 
years, the structure evolution was derived from molecu- 
lar dynamics simulations, where the high computer cost 
imposes the use of sijji™lificd-|Cuptentials based on effec- 
tive medium theory ;Boi3Q'ljElall3 subsequently free elec- 
tron methods were applied to estimate the conduction of 
the metallic neck (or confining potential). More precise 
methods, considering the electronic structure, have also 
been applied to calculate theJNW conductance, but for 
static atomic configurations .til 

Recently, in situ high resolution transmission electron 
microscopy (HRTEM) experiments have provided-a new 
insight in the field. For example, Rodrigues et alx3 have 
showed that just before rupture, gold NW's are crys- 
talline and free of defects and they assume only three 
kinds of atomic arrangements: two of them form bipyra- 
midal constrictions which evolve to one-atom-thick con- 
tacts, while the other one generates rod-like NW's that 
break when they are rather thick|-(jthree-to-four atoms). 
By considering that Onishi et aZJla have already shown 
that atom size contacts display a conductance of 1 Go, 
it is possible to discriminate between the rod-like and 
pyramidal NW morphologies. Further, the relationship 
between each NW type and electrical transport mea- 
surements was obtained by simple crystallographic ar- 
guments. Validation of the initial assumption was ob- 
tained by statistically comparing the occurrence of ob- 
served curve profiles. 

Although these important developments mostly fo- 
cused on the last conductance plateau, a quantitative 
understanding of the correlation between atomic struc- 
ture and conductance during the nanowire stretching is 
still lacking. In this work, we have addressed the con- 
nection between gold NW structure and the quantized 
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conductance behavior during the NW elongation. We 
have used HRTEM to obtain the detailed information 
of the atomic structure evolution of gold NW's during 
stretching, and using crystallographic arguments, we pro- 
posed the three dimensional structure of these nanos- 
tructures. The NW conductance was measured using an 
independent, dedicated experimental set-up: a mechan- 
ically controlled break junction operated in ultra-high- 
vacuum (UHV-MCBJ). In order to correlate both sets 
of data, we have used a semi-empirical atomistic thep=. 
retical technique based on the extended Huckel theorytil 
(EHT) that allows for the estimation of transport prop- 
erties of systems with a large number of atoms J13 The re- 
sults displayed an excellent agreement between observed 
structural and conductance experiments and theoretical 
calculations, enabling a complete modeling of the NW 
elongation process both from the structural and trans- 
port properties points of view. 

The remaining of this work is organized as follows: 
section II describes the experimental techniques used to 
obtain the structural and transport data, respectively; 
section III describes the theoretical model employed to 
calculate the electronic quantum transport through the 
metallic NW's; the experimental results are presented 
in section [j^ and the analysis of the data via theoreti- 
cal analysis is presented in sections [v] for NW's oriented 
along the [100], [111] and [110] directions as well as [v|, 
where the pre sence of impurities is investigated; finally, 
in section 



VII we draw our conclusions. 



II. EXPERIMENTAL APPARATUS 

We have generated NW's in situ in a HRTEM (JEM 
3010 UHR, operated at 300 kV, point resolution 1. 
using the method developed by Takayanagi's group. I 
The procedure consists in focusing the microscope elec- 
tron beam (current density 120 A/cm 2 ) on the sample 
to perforate and grow neighboring holes until a nano- 
metric bridge is formed between two of them. We have 
used a self-supported polycrystalline gold thin film (5 nm 
thick, deposited on a holey carbon grid) in order to gen- 
erate NW's between apexes of differ^Bt orientations and 
elongate them in different directions.!! 2 ! Atomic resolution 
image acquisition has been performed after reducing the 
electron beam intensity to its usual value (30 A/cm 2 ). A 
high sensitivity TV camera (Gatan 622SC, 30 frames/s) 
associated with a conventional video recorder was used 
to register NW real-time evolution. The procedure de- 
scribed above, allows us to generate NW's with a re- 
markable stability, because the NW, its apexes and the 
surrounding thin film regions, form a monolithic block. 
In consequence, the NW's formed by a few atomic layers 
usually show a long life time (1-10 minutes). Despite this 
stability, the generated NW's elongate spontaneously, get 
thinner and then break due to the relative slow move- 
ment of the NW apexes. These apexes displacement are 
probably due to a film deformation induced by thermal 



gradients between parts of the sample, as usually ob- 
served in TEM thin film work. A critical aspect when 
studying such tiny nanostructures is the presence of con- 
taminants, the most critical one in TEM is amorphous 
carbon.EJ However, we must keep in mind, that the in- 
tense electron irradiation transforms carbon into bucky- 
onions and finally clean the gold surface &3 

The electric transport properties of gold NW's were 
studied using an independent instrument special!*-, de- 
signed for this purpose, namely an UHV-MCBJE3 In 
this approach, a macroscopic gold wire (99.99 % pure, 
4> = 75 fj,m) is glued in a flexible substrate in two points: 
then it is weakened in a point between the two fixing 
parts by an incomplete cut. By bending the substrate 
in situ in the UHV, we break the wire and produce two 
clean gold surfaces; using the same bending movement, 
the fresh tips are put together and separated repeatedly 
in order to generated and deform NW's. It's important 
to remark that in this configuration the NW's are gen- 
erated from surfaces obtained in UHV (< 10~ 8 Pa), pa 
it is expected to have a clean sample for few hours .1121 
This care with vacuum conditions is extremely impor- 
tant to generate-jjeliable experimental measurements on 
NW properties J3EI 

Electrical measurements were done using a home-made 
voltage source, and current-voltage converter powered 
with isolated batteries to reduce electrical noise. Data 
acquisition was based on an eight-bit digital oscillo- 
scope (Tektronic TDS540C, 2 GSample/s, bandwidth 
500 MHz); this electronic system was developed with the 
aim of improving the conductance measurement preci- 
sion, yielding values with a relative error of (AG/G) ~ 
10~ 4 . The applied voltage was 100 mV, and the con- 
ductance was measured in the [0, 2.8] Go range to im- 
prove the detection of the last two quantum conductance 
plateaus. 



III. THEORY 

We follow the approach introduced by Emberly and 
Kirczenowli3 that is based on the extended Huckel the- 
ory, which is a method traditionally used in quantum 
chemistry to calcuLate-the electronic properties of clus- 
ters and molecules. tSEil One advantage of this approach 
lies in the possibility of treating the whole nanocontact - 
both tips and the nanowire between them - as a unique 
entity, described by molecular orbitals. 

The tips and the nanowire, upon which we base our 
transport calculations, are built on the basis of the crys- 
tallographic properties of the material as well as the ob- 
tained HRTEM images; no attempt was made to sim- 
ulate the dynamics of the NW formation. In our cal- 
culations the nanocontacts are made up of as many as 
150 gold atoms for some structures. In order to calcu- 
late the quantum electronic states of Au clusters as big 
as these, a semi-empirical method is necessary to pro- 
duce results within reasonable time. For this reason the 
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EHT is employed; the method includes overlap terms and 
matrix elements that are not restricted to first neigh- 
bors; moreover, all the valence states - 5d, 6s and 6p, 
in the case of Au - are considered. One drawback is 
the fact that, in its standard form, the EHT is a one- 
electron formalism that does not take into account ex- 
plicitly electron-electron interaction effects. However, we 
believe that this problem can be partially circumvented, 
because of the parameterization of the orbital wavefunc- 
tions by self-consistent Hartree-Fock states and the use 
of semi-empirical energy parameters. Such effects are 
more relevant for the analysis of the structural stability 
and relaxation of nanowires, where the calculation of en- 
ergy manifolds is necessary. Nonetheless, tight-binding 
together with molecular dynamics methods, similar to 
the present one, have been successfully used to simulate 
the dynamics of the nanowire formationE3 



cells whose MO's are also obtained by the EHT method: 



Left Lead 



Right Lea d 




FIG. 1: Sketch of model structure used in the calculations 
evidencing one unit-cell of the incoming and outgoing leads 
and the actual NW. 



Let us assume that a nanocontact is formed by two 
apexes (or tips) and the nanowire itself that connects 
them. We start the description by defining the molecular 
Hamiltonian of the nanocontact, Hnc and its molecular 
orbitals (MO's) |</>): 
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(1) 

(2) 



In (0) the sum over v extends over the 5d, 6s and 6p 
orbitals and the index j covers all the Au atoms in the 
nanocontact. The molecular orbitals \<f>), written as a 
linear combination of atomic orbitals in (hi), obey the 



orthonormality relation 



Equation (H) de- 



fines the MO energies. So that these MO's bear quantum 
phase coherence along the entire nanocontact. 

To account for electron transport the nanocontact is 
coupled to infinite Au leads on both sides (left (£) and 
right (TZ) of the nanocontact), which behave as charge 
reservoirs (see Fig. |l|). The leads are divided into unit 
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(3) 
(4) 



For the sake of simplicity only 6s Au orbitals are used 
in (|^). The Hamiltonian H^r-jz) includes overlap terms 
and is not restricted to first neighbor matrix elements 
within the unit cell. We consider this approximation be- 
cause the contacts are expected to behave as bulk gold, 
hence s orbitals are adequate for describing its electronic 
properties. In the nanoconstriction, where the number 
of atoms is smaller, we expect a different behavior. Al- 
though (a'\a) = 5 a i a , it is noteworthy that the scalar 
product between both types of MO's define the overlap 
term: 
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(5) 



Having defined the basic physical quantities, let us con- 
sider the total Hamiltonian 



H = H C + Hue +H n + Wc + W n ■ 



(6) 



The term Wc in (Q) is responsible for coupling the 
nanocontact with the first unit cell of the lead on its 
left and W-jz, similarly, with respect to the first unit cell 
on the right. In principle He and H-r can be different 
from each other, as in the case of tips with distinct crys- 
tallographic orientations or defects. 

To account for the transport calculations consider 
Bloch electrons that propagate through the leads: 



l*> = 



E< 



|/3 n) 



with |/3 n) as one of the MO's of the unit cell at position 
n and |/3) G {|a)}- At the site of the nanocontact the 
incoming Bloch wave is scattered giving rise to reflected 
and transmitted waves, according to ihe framework of 
the multi-channel Landauer formalism. E3 Therefore, for a 
rightward incoming electron we write the following wave- 
function |*) = + \^mc) + |*72.), where 
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(10) 



The incoming electron has energy E(0p) > in order 
to be a propagating state. In expression @, \<j> 0) rep- 
resents the </)-MO's of the nanocontact, which occupies 
the site n = 0, and C^,p are unknown coefficients. Sim- 
ilarly, in M) and (|To|), \a n) is the a-MO of the unit 
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cell at position n\ r Q/ 3 and t a p are the channel reflection 
and transmission coefficients, to be determined in the se- 
quence, which will be used to calculate the total transmis- 
sion T(E) and reflection R(E) of electrons propagating 
through the nanocontact. The parameter 8 is a dimcn- 
sionless wave- vector for the Bloch waves in the leads. 

In order to determine the coefficients {r a f3,t a /3,C c j >a } 
a complete set of MO's {£} = {(f)} U {a} is put together, 
where we project the equation H\^f) — E\^). As a result 
we obtain a non-homogeneous system of coupled equa- 
tions that can be solved numerically to yield these co- 
efficients. A detailed description is provided in the ap- 
pendix. 

Concluding the procedure, we define the total trans- 
mission and reflection coefficients, which are written in 
terms of the incoming and outgoing coefficients and the 
electron velocity v a of the respective Bloch wave mode: 



R,(E) = J2 ~M 



(ii) 

(12) 



where the prime means that the sum is made exclusively 
over the propagating modes, tha44s, E(9 a ) > (a- 

Finally, the Landauer formulaofiJ is used to determine 
the charge transport through the nanocontact. For a 
degenerate metallic system in the linear response regime, 
the electrochemical potential difference between source 
and drain takes the limit A[i — > and the conductance 
G = I/V becomes 



(13) 



where r\ is the Fermi-Dirac distribution and Ep is the 
Fermi energy. Since the voltage applied across the leads 
is small compared to Ep the conductance is to be con- 
sidered when the electron energy is close to this value. 



IV. EXPERIMENTAL RESULTS 

Figure ||(a) shows the global histogram generated from 
a series of 500 collected conductance curves of gold NW's 
measured using the UHV-MCBJ. Three peaks can be 
identified in the histogram; as expected the major one 
is close to 1 Go and, the two minor ones are at 1.8 and 
2.8 Go- r-i 

Rodrigues et al£3 have demonstrated that analyzing 
the occurrence of the last two conductance plateaus 
(around 1 and 2 Go), the conductance curves can be 
grouped in three different classes: a) only a plateau at 1 
Go b) plateaus at 1 and 2 Go; and, finally, c) no plateau at 
1 Go- These characteristics are associated with the crys- 
tallographic direction along which the NW is being elon- 
gated, corresponding to the [100], [111] and [110], respec- 
tively. As the global histogram superposes and hides this 
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FIG. 2: (a) Au NW global histogram and (b-d) histograms 
associated with each of the three types of curve profile (see 
text for explanation). The greatest recorded value shown in 
the histogram's abscissa is 2.8 Go- The arrows show the ex- 
istence of peaks below 3 Go- The insets of every figure show 
typical conductance curves for each direction. 



structural information, we have built three histograms 
formed by each group of conductance curves (Fig. |^(b- 
d)), examples of the conductance curves are presented at 
the insets). This approach allows us to obtain a statis- 
tical validation for the sequence of conductance plateaus 
appearing during the NW elongation along each axis, ex- 
tending the range from 2 to 3 Go- In these terms, it can 
be deduced that the thinning of a NW along [100] axis 
(hereafter denoted [100] NW) generates a conductance 
curve with plateaus close to 3 and 1 Go (see curves inset 
Fig. 2(b)); a [111] NW, plateaus around 3, 2 and 1 G 
and; finally, for [110] NW's, plateaus around 3 and 2 Gq. 




FIG. 3: Real time HRETM images showing typical nanowires 
formed in the direction of the three zone axes, (a-d) Elonga- 
tion process of a [100]-oriented NW; note the existence of two 
misaligned planes just before the formation of the ATC. (e) 
ATC oriented along the [111] direction and (f) rod-like struc- 
ture formed along the [110] direction. 



5 



Figure || shows typical HRTEM micrographs of the 
three possible Au NW structures; a common feature of 
all the observed NW's is the fact that they are defect free 
crystals. Figures ||(a-d) show snapshots of the elongation 
process of a NW along the [100] direction, where the pro- 
gressive thinning of the constriction can be followed with 
atomic resolution. We can see from these images that the 
NW stretches from a symmetrical bipyramidal NW with 
a two- atom wide minimal section (Fig. |](b)) to a con- 
striction formed by two misaligned truncated pyramids 
with a two atom wide top (Fig. ||(c)). Finally a one- 
atom-thick contact is formed (Fig. f|(d)); this contact 
subsequently evolves to form a suspended chain of atoms 
(ATC) before rupture (not shown). Figure ||(e) allows 
the clear visualization of the bipyramidal morphology of 
a one-atom-thick contact formed for [111] NW's. In con- 
trast, [HOlNW's generate rod-like structures before rup- 
ture (Fig. |(f)). 

At this point, it is useful to remind some basic aspects 
of the reported in situ HRTEM experiments. Firstly, 
HRTEM images of NW's can be described as a two- 
dimensional projection of the atomic potential, then good 
quality atomic resolution information is only obtained 
from crystals very well oriented in relation to the micro- 
scope electron beam. In our experiments, we can not 
easily control the orientation of the apexes forming the 
nanometric constrictions, so the possibility of getting se- 
quences as shown in Fig. |(a-d) ([100] NW) is statistical 
and unpredictable. 



V. DISCUSSION: CORRELATION BETWEEN 
ATOMIC STRUCTURE AND ELECTRONIC 
EFFECTS 

A. [100] Gold NW's 



a) 941 ( [100] Tip ) 



! b) 949 
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FIG. 4: Atomic geometry of different Au nanocontacts ori- 
ented along the [100] direction. The notation indicates the 
number of atoms in the planes perpendicular to the [100] axis 
in the constriction region: (a) 941, (b) 949 and (c) 9449 tips. 
The dashed lines in (a) indicate the edge of the fee unit-cell 
and ao is the lattice parameter. 



Unfortunately, this kind of detailed information on the 
atomic arrangement evolution could not be observed for 
the other two NW types of nanowires ([111] and [110]). 
Another limiting factor of our dynamical HRTEM obser- 
vations is the temporal sampling of our image acquisition 
system allowing the capture of 30 frames per second; very 
rapid events (< 1/30 s) may go by undetected. 

Despite these difficulties, HRTEM images provide a 
wealth of remarkable information to derive the atomic 
structure of gold NW's from. By combining the atomic 
resolution micrographs (two dimensional projections of 
the structure) and the known surface energy properties 
of the analyzed systems it is possible to deduce the three 
dimensional atomic arrangement of NW's. This proce- 
dure is based on the so-called Wulff construction that 
allows a quick geometricalrjdetermination of the faceting 
pattern of nanostructures;CJ such approach .has already, 
been successfully applied to study Aujij AgtS and PtEZl 
NW's. 

In the following sections, we present a careful and thor- 
ough analysis of the derived atomic arrangements and 
their evolution for the three kinds of Au NW's; the theo- 
retically obtained electronic transport properties of such 
structures are compared with the conductance behavior 
observed with the UHV-MCBJ experiment. 



The [100] nanocontacts are generated by two pyrami- 
dal tips, each of them faceted by four low energy (111) 
planes; they can be elongated until atomic-chains a*s 
formed, having generally two-to-four atoms in length.tHI 
Figure ^ illustrates three nanocontact geometries based 
on the crystallographic structure of face-centered-cubic 
(fee) bulk gold and the Wulff method. These structures 
simulate the steps of the NW evolution in the exper- 
iments. Figure [l](a) shows the cross-section ([100] di- 
rection) and side view ([011] direction) of an atomically 
sharp [100] NW tip; the side view presented represents 
the structural arrangement as seen from HRTEM images 
in Fig. ||(a-d). This particular tip is formed by the 
stacking of a plane of 9 atoms (black circles) followed 
by a plane of 4 atoms (dark gray circles) and, finally, one 
atom at the end of the tip (light gray circle). The nota- 
tion used hereafter to designate the structures is given by 
the number of atoms that make up the alternate planes 
stacked along the NW orientation axis. The one-atom- 
thick contact (941149) which is the structure observed 
just before rupture (example is presented in Fig. |3|(e)) 
is then formed by mirror imaging the tip arrangement 
described above. We can use this pyramidal tip arrange- 
ment to model the thicker NW's observed prior to the 
formation of the 941149 structure, i. e., we can build 
thicker structures by removing planes of atoms from the 
tip arrangement in Fig. 0(a). Hence, in simple terms, we 



6 



would expect that the second thinnest [100] constriction 
would be formed by pyramids where the one atom tip 
is removed leading to a minimum cross-section of four 
atoms. By joining two pyramids sharing the four-atom 
plane we generate the 949 neck shown in Fig. ^(b) and 
corresponding to the experimental results in Fig. [§(b). 

However, HRTEM data clearly show the existence of 
an intermediate structure (Fig. |3|(c)) where the two trun- 
cated pyramids forming the apexes are not aligned with 
respect to each other. Assuming that all the neck region 
is defect free (i. e. a perfect arrangement of atoms in a 
fee lattice) and following the Wullf morphology we must 
conclude that the NW is formed by two four-atom planes 
which are not not aligned with each other, rather the cor- 
ner of one of these planes is located at the central axis of 
the other and vice-versa. This structure is presented in 
Fig. ||(c) and labeled 9449. 
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FIG. 5: Calculations of the quantum conductance (G) as a 
function of the electron energy for various NW structures ori- 
ented along the [100] direction: 941149 (the one- atom-thick 
contact), the 9449 and 949 nanocontacts. The vertical line 
indicates the Fermi energy; the conductance is in units of the 
quantum Go. 

Using the formalism of section III, we have calcu- 
lated the quantum conductance for the NW morpholo- 
gies shown in Fig. |as a function of the electron energy. 
The results are shown in Fig. |^, where the vertical line 
indicates the value of the Fermi energy in the leads. At 
first, we notice that the one-atom-thick contact (941149, 
distance between the atoms is 2.88 A even at the con- 
tact) shows a conductance curve that yields a very sta- 
ble plateau at G(Ep) = 1 Go which corresponds to a 
peak around 1 Go in the histograms (Fig. ||(b)). On the 
other hand, the 9449 nanocontact display a conductance 
of G{E F ) ~ 3 G , and the 949, 4.5 G . 

In brief, the one-atom-thick (941149) and the ATC's 
conductance (9411149, also calculated but not shown 
here) is ~ 1 Go while the 9449 structure shows a con- 
ductance of ~ 3 Go- Each of these conductance values 
agrees very well with the peaks in the global conduc- 
tance histogram for the [100] zone axis (Fig. |2|(b)). In 
other words, he have accounted for the peaks in the his- 
togram corroborating the assumption that the sequence 
of atomic arrangements shown in Fig. |4| can, in fact, 



represent the dynamics of NW evolution along the [100] 
direction. It is also worth noticing that the deduced 
structural evolution shows no contribution around 2 Go, 
in agreement with the conductance histogram. Finally, 
the 949 morphology presents G(E F ) ~ 4.5 Go, which 
would represent a thicker structure however, our exper- 
iments have not included conductance measurements in 
that range. 



B. [Ill] Gold NW's 
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FIG. 6: Atomic geometry of different Au nanocontacts ori- 
ented along the [111] direction. Figure (a) shows the arrange- 
ment of the tip. The notation indicates the number of atoms 
in the planes perpendicular to the [111] axis in the constric- 
tion region: (b) 731137, (c) 737, (d) 7337, and (e) 732237 

For the case of [111] NW's, we have no atomic resolu- 
tion HRTEM images that show the structural evolution 
during the thinning process. However, we can modc i-tlic 
structures using the same procedure of section V_A,li3 i. 
e., we can start with the one-atom-thick nanowires and 
make our way up towards thicker structures. Again, the 
tip morphology was built from a fee crystal. It is worth 
noting that the stacking of planes in the [111] direction is 
made by three hexagonal planes placed at alternate inter- 
site spaces (usually labeled ABCABC. . .). The derived 
[111] pyramidal morphology is presented in Fig. ||(a). 
This arrangement is obtained from the smallest possi- 
ble fee gold nanoparticle, i. e., the cuboctahedral pith 
regular hexagonal faces based on the Wulff methodO 

As it was done in the previous section, the one-atom- 
thick NW is formed by two aligned [111] apexes oriented 
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in opposite directions. We can label this structure 731137 
by following the notation introduced earlier. This struc- 
ture is depicted in Fig. ||(b). In analogy with the pro- 
cedure adopted for the [100] direction, we would initially 
expect that the removal of the last atom from the tip and 
then by joining the remaining truncated pyramids we 
would have the arrangement before the one-atom-thick 
NW is formed. This particular arrangement is sketched 
in Fig. ^(c) and is named 737. However, this arrange- 
ment is defective because it forms a twin at the neck. 
This twin should be unfavorable in the|-sce*c-,that the 
NW's are free of defects in this size range JI3£§El. In that 
case, we proposed a structure analo gous to the 9449 [100] 
NW which was presented in section |V A| . This new struc- 
ture is formed by adjoining two truncated [111] pyramids 
(73) rotated 180° with respect to each other. The result- 
ing arrangement has a perfect fee stacking as shown in 
Fig. ^(d). Finally, molecular dynamics simulations using 
classical potentialsEHl indicate the existence of a thinner 
region with two misaligned dimmers. As a matter of fact, 
this arrangement is formed by simply removing atoms 
from the apexes of the 7337 NW. This structure, labeled 
732237, is presented in Fig. ||(e) for atoms in a perfect 
fee arrangement. 



in Fig. [7| Conductance measurements for this orienta- 
tion (Fig. 0(c)), show the existence of a peak close to 2 
G - In this sense, the 732237 NW, results in G{E F ) ~ 
1.9 Go as can be seen in Fig. (?]. Finally, structure 731137 
yields G = 1 Go (Fig. 0(d)). 

One might be led to believe that the relaxation of the 
thinnest region (..22..) of the 732237 structure will ro- 
tate the two dimmers to form a perfect tetrahedron in 
order to maximize the number of neighbors. This defor- 
mation generates a NW with two [111] apexes bounding a 
[100]-oriented neck. This deformation decreases the con- 
ductance by a significant amount from 1.9 to ~ 1.2 Go; 
emphasizing that nanowires are particularly sensitive to 
breaking symmetries and to defects in the neck region. 

In other words, the 7337 morphology has a conduction 
close to 2.6 G , the 732237 NW shows the calculated 
conductance value G{Ep) ~ 1.9 Go and at last, for the 
one-atom-thick contact, G(Ep) ~ 1 Go- Each of these 
structures can be associated with one of the peaks of the 
global conductance histogram of Fig. §(c), and to the 
evolution steps of nanowires along this direction. 



C. [110] Gold NW's 
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FIG. 7: Calculations of the quantum conductance G as a 
function of the electron energy for various proposed nanocon- 
tact structures oriented along the [111] direction: the 737, 
7337, 732237 and the one-atom-thick (731137) nanocontacts. 

In Fig. [7|, we present the theoretical conductance cal- 
culations for each of the proposed structures for [111] 
NW's. As previously stated, the experimental results 
must be compared to the calculated conductance at the 
Fermi energy level, represented by the vertical line in 
the figure. Although the 737 arrangement (Fig. ||(c)) 
would be a possible configuration displaying a conduc- 
tance of ~ 2.8 Go, we must bear in mind that it presents 
a twin defect. Hence we have neglected it in our analysis. 
The calculated conductance for the 7337 structure, which 
presents a perfect fee stacking, is below 2.6 Go as shown 
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FIG. 8: Cross-section (left) and side view (right) of the 
schematic representation of possible rod- like [110] NW's at 
various stages of the thinning process. The atomic arrange- 
ments are labeled from the number of atoms forming the alter- 
nate stacking planes generating the rods: (a) 4/3, (b) 2/2, (c) 
1/1 and (d) 2/0. See the text for a definition of the notation. 

The HRTEM results show that compared to the other 
two kinds of nanowires, the [110] NW's are rather brittle, 
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i. e., they break when they are two-to-three atoms thick 
without forming atomic contacts. Figures ||(a-d) show 
the proposed structure of rod-like [110] NWis, as they 
decrease in diameter, by loosing atomic planes&j. Again, 
the notation used in naming these structures comes from 
the number of atoms contained in each plane forming 
the NW (black and gray circles). However, in this case, 
due to their large aspect ratio, we show only the rod- 
like part and leave out the apexes. The thickest NW 
(Fig. ||(a)), labeled 4/3, is formed by two families of low- 
energy (111) facets and two high-energy (100) ones. By 
removing 3 lines of atoms from the previous structure 
we reduce it to the 2/2 NW. Finally we consider two 
possible scenarios: the 2/2 NW is subtracted of one of 
the (111) planes (atoms Ai and A2), giving rise to the 
alternate 2-atom chain, 1/1, shown in Fig. ||(c)); the 
other possibility is the elimination of the central plane of 
atoms (atoms A2 and A4), generating the parallel 2-atom 
chain 2/0 (Fig. |(d)). 



We might consider that the elongation process will in- 
crease the interplanar distance in the 2/2 structure and 
bring atoms A2 and A4 together until they reach first 
neighbor distances. Nonetheless the orientation of the 
stacking is changed to the [100] dir ection . As it happened 
for the 732237 structure of section V B , the conductance 
calculated for this case shows slight decrease compared 
to the perfect structure (from 2 to 1.5 Go). 

In the light of the discussion above, we can model the 
atomic rearrangements in this direction by removing lines 
of atoms, starting at the 4/3 structure, which has a con- 
ductance G{Ep) > 4.5 Go, down to the 2/2 that can be 
associated, aided by the calculations, to the peak at 2.8 
Go on Fig. |](d). Finally, the peak at 2.0 Go is related 
to either the 1/1 or 2/0 structures. 



VI. EFFECT OF CARBON IMPURITIES 
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FIG. 9: Calculations of the quantum conductance G as a func- 
tion of the electron energy for various rod-like NW's along the 
[110] direction: (a) 4/3, (b) 2/2, (c) 1/1, (d) 2/0 structures. 
The vertical line indicates the Ef- 



The calculated conductance for each of the [110] NWs 
presented in Fig. || is shown in Fig. ^. We can see that 
the conductance curves presented show oscillations that 
arise from interference effects inside the contacts. This 
oscillations are sensitive to the atomic positions so we 
must consider an average around the value E^. The cal- 
culated conductance for the 4/3 structure, as presented 
in Fig. |(a), is G{E F ) > 4.5 G . Meanwhile, for the 2/2 
case it is close to 2.8 Go- The last two NW's structures, 
1/1 and 2/0, show essentially the same conductance pat- 
tern, G(E F ) ~ 2 G . The last HRTEM image - within the 
time resolution considered - for NWs formed along the 
[110] direction is the 2/2 structure, so it is not possible 
to resolve the final structure. However, the conductance 
histogram of Fig. |2|(d) shows that the last plateau is ap- 
proximately 1.8 Go, indicating that at least one of these 
structures must be present in UHV-MCBJ experiments. 




0.2 1 2 

Conductance [ G ] 

FIG. 10: Conductance histogram of a set of 500 curves mea- 
sured using the UHV-MCBJ after several hours of experiment. 
The arrow shows the appearance of a peak between 0.1 and 0.2 
Go which is not present in the first hours of experiment. Inset: 
Typical conductance curves; note the existence of plateaus at 
G ~ 0.2 Go 



The statistical correlation between the occurrence of 
conductance curve profiles and the multiplicity of zone 
axis has been hitherto described as a property of the Au 
nanowires.|__l However, the accordance is only observed 
during the first hours of measurement after the macro- 
scopic wire was fractured inside the MCBJ, while the 
gold surfaces are free of impurities; passed this initial 
period, contamination occurs even in UHV conditions. 
That is then evidenced by the appearance of conduc- 
tance plateaus at values ranging approximately from QJ. 
to 0.2 Go and whose occurrence increases with time.Ej 
Similar observations have been reported fox. NWs seen 
using UHV scanning tunneling microscopy.t__ The global 
histogram and typical experimental conductance curves 
illustrating this effect are shown in Fig. |l^. Likewise, an- 
other indication that the MCBJ might be contaminated 
is the fact that the 1 Go plateaus become exceedingly sta- 
ble and long, suggesting the formation of entire nanowires 
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made of impurity atoms. One of the most .probable con- 
taminants is carbon. Ab initio calculationstj have shown 
that the carbon atom naturally occupies a position be- 
tween two gold atoms in ATC's. Since the presence of 
carbon cannot be detected-by HRTEM images because 
of its low atomic numberJ^J the effect characterized by 
plateaus at 0.1-0.2 Go (Fig. |l0| ) is therefore the only 
experimental evidence of its presence. 
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FIG. 11: Theoretical results for the conductance of ATC's 
as a function of the energy. The solid line corresponds to a 
pure gold ATC (interatomic distance 2.88A). The dashed line 
shows the effect of carbon impurities in the ATC, as depicted 
in the figure. The Au-C interatomic dptance is fixed at 1.8 
A, obtained from ab initio calculations E3 

We have used the formalism of section III to address 
this question. For this purpose, consider the following 
gold ATC's in the [100] direction: a chain of 2 suspended 
gold atoms between two [100] tips (Au-Au-Au-Au, four 
Au atoms forming the narrowest region), as well as the 
contaminated Au-C-Au-Au chain. The, assumed distance 
separating C and Au atoms is 1.8AJ23 whereas the Au- 
Au distance is taken to be 2.88A. Calculations of the 
conductance produced for the 2 cases are illustrated in 
Fig. [H]. It can be verified that the presence of C atoms 
in the chain is responsible for decreasing the conductance 
of the nanocontact considerably from 1 to ~ 0.1 Go, in 
very good agreement with the experimental results and, 
moreover, corroborating the assumption that carbon can 
act as a contaminant of Au NW's after a few hours of 
measurement. 



VII. CONCLUSION 

In summary, we have studied experimentally and the- 
oretically the elongation of gold NW's, both from the 
structural evolution and the quantum conductance be- 
havior. From the experimental point of view, we have 



used two independent set-ups: an UHV-MCBJ operat- 
ing at room temperature for conductance measurements 
and HRTEM for real-time imaging of the NW atomic 
rearrangements. 

Gold NW's adopt only three kinds of atomic struc- 
tures where the elongation direction is parallel to one of 
the [100], [111] and [110] axis. For each one of these NW 
types, we have derived the three dimensional structures 
during different stages of the elongation process, by com- 
bining HRTEM images and the Wulff method. Subse- 
quently conductance calculations, based on the extended 
Hiickel theory and the Landauer formalism, were applied 
to the deduced atomic structures, showing an excellent 
accordance with the conductance measurements. This 
approach has provided a consistent correlation between 
the structural evolution of the three kinds of Au NW's 
and the conductance behavior in the 0-3 Go conductance 
range. 

The EHT, although rather simple, allows us to take 
into account the atomic arrangement of the nanocontacs 
and, also it includes the possibility of considering atoms 
of different species. As an example of its versatility, we 
have considered carbon impurities on suspended atom 
chains; the results indicate that they act as a contami- 
nant causing the conductance plateaus to decrease from 
1 to 0.1 Go- This result agrees extremely well with our 
conductance data of coptaminated Au NW's and, also 
those from Olesen et alx3 

We have, therefore, thoroughly analyzed the relation- 
ship between the structural modifications involved in the 
elongation of Au nanowires and their transport behav- 
ior. The main conclusion that must be drawn from this 
work is the determinant role of an atomistic description 
when analyzing the conductance properties of nanoscopic 
metallic wires. 
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IX. APPENDIX 

This appendix presents the details of the theory, cir- 
cumvented in section III. When projected onto the basis 
states (£V| the equation H\%) = E\^) yields the fol- 
lowing set of coupled equations for the set of coefficients 

V ce/3) ta,P> G<jj j( j}: 



J 



E E + E + E E = - E a 4 ) 

n<0 a (f> n>0 a n<0 
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with 



4'4 = (t'n'\H - E\£n) 



(15) 
(16) 



The term H^A renders the MO energies of a given unit 
cell whenever £' = £ and n' — n. It also accounts 
for the hoppingenergies between the nanocontact and 
the first unit cell of the leads, on both sides, that is, 
(<j>'n'\H\ari) = Wqj°[. The same idea applies to the 

overlap terms 55/1 • Within the leads we simplify the 
treatment by taking into account only the coupling be- 
tween MO's of the same type, which yields (a'n'\H\an) = 
h 5 a ' a 5 n ' ,n±i and (a'n'\an) = Sq 5 a 'a&ri ,n±i, with h and 
So being the same for all MO's. This is necessary for ob- 



taining analytical expressions for the dispersion energies 
of the electrons in the leads, which allows us to consider 
explicitly only the sites n = —1,0,1 in (|l4|). Thus the 
energy dispersion relation is 



E{9 a ) = 



e a + 2h cos (8 a ) 
1 + 2S cos (0 a ) 



(17) 



and the velocity of the electron mode v a = dE a /dk ai 
with k a = 9 a /a and a is the length of the unit cell. The 
parameters h = 4 eV and Sq — 0.15 were chosen to yield 
a bandwidth of approximately 8 eV, which is found in 
bulk Au. 

Finally we obtain for the representative sites : 



n' = -1 



n' = 



ri = 1 



Y / {(h-ES )e i ^ + (ei-E)e iS "}S^' r aff +J2(W^t - ES^*)C^ 

Q (f) 

= -{(h- ES )e- m ? + (e% - E)e-^}5% , 



(18) 



a (pa 

^2 {(h - ES Q y 2e « + & - E)e^}8i tap + ^(Wfo 4, - ESf*)C^ = 



(19) 



(20) 
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